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A route to quasi-perfect invisibility cylindrical cloaks without extreme values in the 

parameters 

Shcng Xi, Hongshcng Chcnp Bailc Zhang. Bac-Ian Wu, and Jin Au Kong 
Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA, 
and The Electromagnetics Academy at Zhejiang University, Zhejiang University, Hangzhou 310027, China 

The method of coordinate transformation offers a way to realize perfect cloaks, but provides less 
ability to characterize the performance of a multilayered cloak in practice. Here, we propose an 
analytical model to predict the performance of a multilayered cylindrical cloak, based on which, the 
cloak in practice can be optimized to diminish the intrinsic scatterings caused by discretization and 
simplification. Extremely low scattering or "quasi-perfect invisibility" can be achieved with only a 
few layers of anisotropic metamaterials without following the transformation method. Meanwhile, 
the permittivity and permeability parameters of the layers are relatively small, which is a remarkable 
advantage of our approach. 



PACS numbers: 41.20.Jb, 42.25.Fx 
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Various efforts have been made on the reahzation of invisibility [1|, |2, |3|, U, |5[. Pendry et al. theoretically proposed 
the perfect invisibility cloak for electromagnetic waves l| , utilizing anisotropic and inhomogeneous media to mimic the 
space squeezing. Later, the effectiveness of the transformation based cloak was demonstrated by ray tracing[6j, full 



IC 



11 



12| , as well. In practice, the difficulty in 



wave finite element simulations |7|,l8|, and analytical scattering modelsj; 
construction a perfect invisibility cylindrical cloak is the requirement of continuous inhomogeneity and high anisotropy 
with extreme values in the parameters. Simplified parameters based on the coordinate transformation were then 

r ~ 

utilized to facilitate the physical realization 



iimplified 

iBQ,[i4i 



15j. Constraints 



in expense of the aroused inherent scatterings . 
on the bandwidth were studied as well[16j. The first sample of cylindrical cloak has been created usin g m ultilayered 
metamaterials 



17j . but a lot 



13[ . Bi-layered isotropic media was also proposed for achieving the effective anisotropy . 
of thin layers are needed which increases the construction complexity. Moreover, the transformation method provide 
less ability to predict the performance of a practical construction composed of discontinuous layers of homogeneous 



anisotropic metamaterials. Therefore, it is very necessary to investigate a better way to design a practical cloak with 



'Electronic address; lchenhs@ewt.mit.edul 



2 




FIG. 1: Configuration of a multilayered cylindrical cloak. 

good performance. 

In this paper, in order to get the exact behavior of a multilayered cloak, the analytical model of a cylindrical cloak 
created with multilayered anisotropic materials is established based on the full wave scattering theory 
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igf. Oi 



results show that, by using only a few layers of anisotropic materials, a "quasi-perfect invisibility" multilayered cloak 
with near zero scattering can still be achieved without following the design method of coordinate transformation, 
and the parameters obtained are relatively small and possible to be realized by metamaterials. The impedances 
between the adjoined layers do not really match each other but can produce zero reflection, which can be treated as 
the counterpart in cylindrical geometry of the reflcctionlcss one-dimensional multilayered slab. All of these results 
provide a second and better way of designing a multilayered cloak. 

We use cylindrical cloak as an example. Without losing the generality, the case a TE-polarized plane wave with 
unit magnitude normally incident onto an A/-layer cylindrical cloak (from region 1 to M) is considered, as shown 
in Fig. [1] The TM case can be analyzed similarly. The radiuses of the boundaries of the cloak arc denoted by Rm 
(m = 0, 1, • • • , M). The relative constitutive parameters in region m arc assumed to be constants denoted by fipm, 
fi^rn and ezm while the region m = is assumed to be free space and the core region {p < Rm) is assumed to be 
PEC. The electric fields E^m in the region m satisfy the following equation, 



3 



1 d f p dE^,n\ I d f I dE, 



pdp\p.^„i dp J p'^d(j)\ppm d(j> 



''^O^zmEzm — 0- (1) 



By applying the method of separation of variables, the general expression for the electric fields in region m can be 
expressed as 

oo 

Ezra — / ^ ^van y-^VmnO^i^P) ^m{m+l)n-^i^mn 

exp(m(/)), (2) 



n— — oo 



where Vmn = ny^ p^m/ Ppm and the wave number in region m is km = fco^ezmM^m- Different from the isotropic 
layered case, here, Vmn is a fraction. Jum„ > ^i^mn represent the i^m„th order Besscl functions of the first kind and 
the t^m„th order Hankel functions of the first kind, respectively. a,„„ is the unknown coefRcients and r^(,„+i)„ 
is the scattering coefficient on the boundary R„i- When a standing wave incident from region m onto the the 
boundary the direct reflection coefficient, which represents the ratio between the directly reflected wave and 
the incident wave, is r„i{m+i)n — U'.h ~ Vm/ilm+ijji)/ i^h' ji + rjm I ilm+ihji) , and the direct transmission coeffi- 
cient, which represents the ratio between the directly transmitted wave in region jti + 1 and the incident wave in 
region to, is tm(m+i)n = — 2z/(7rfc„ii?,„)/(— /I'ji + rjm / rjm+ihi[) . Similarly, for an outgoing wave the direct reflection 
and transmission coefficients on i?„ are r(^m+i)mn = [K^ - ilm+i/rimhih')/{-j[h + ri,n+i/ilmjih') and = 

2i/{Trkm+lRm)/i~j'lh + r]„i+l/'nrnjlh'). Here j = Jv^^ikmRm), j' = jl^^^ikmRm), jl = Jv^^ + i)„{Kn+lRm), 
j'l = J,,^ ^ ^ -^.j (km+lRm)i h = H-u^^^{k„iRm), h' = H^^^^{k„iRm) , ^1 = -ff-U(„,_,_i)„ (^m+l -Rm ) , h[ — H,,^ ^ (km+lRm), 



and 77,„ = yjp^^l^ 



19|. Therefore, the scattering coefficient in layer m (m=0,l,...,M-l) can be written as 



19| 



i(m+l)n ^m(m+l)n ~i~ ^(m+l)mn; (3) 



where = im(m+i)ni(m+i)mn^'(m+i)(m+2)n/(l " »'(m+i)m,i^(m+i)(rn+2)n) represents the wave coming out 

from Rm due to the multiple reflections and transmissions on the boundaries inside i?,„. At fA/(M+i)n = 

~ JvMn^MRu) I iivM-n.^MRM)^ thcreforc all the fm(m+i)n can be derived using Eq. [Sj and a„m can also be derived 
by matching the boundary conditions. The coefficients of the scattering flelds in region are 6o„ ~ aQnXmn- Based 
on the cylindrical scattering model, the far-field total scattering efficiency or the scattering cross section normalized 
by the geometrical cross section for the multilayered cylindrical cloak is obtained as 
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FIG. 2: Comparison of the differential scattering efficiencies when a TE-polarized plane wave normally incident on to the PEC 
cylinder without cloak (dash dotted line) and with the cloak proposed in Ref. [13|] (dotted line). 

In practice, the ideal parameters obtained from the transformation method need to be discretized for realization, 
which will destroy the perfect invisibility of the cloak. Using the proposed method, such effects of descretization 



13| 



and simplification of the transformation based (TB) cloak can be quantitatively analyzed. For example. Ref. 
proposed a 10-layer simplified cloak with /i^ = 1 for the experiment, in which the copper cylinder core with radius 
0.709A is coated by a multilayered cloak with inner radius 0.768A and outer radius 1.670A. Utilizing the parameters 
including the losses in the metamaterials in Ref. 



13|, and assuming the core to be PEC for convenience, for the 
normal incidence of a TE-polarized plane wave, the forward and backward scatterings have been reduced by about 
4.8dB and 4.1dB, respectively, comparing with the bared PEC core, as shown by the dotted line in Fig. [S] The near 
field distributions can also be calculated using our method. Fig. [3] shows the electric field distributions of the above 
case. Our analytical model shows some qualitative agreement with the experimental field distributions shown in 13l |. 



where both reduced forward and backward scatterings can be observed. 

A more interesting thing is that by using the genetic algorithms, a widely used optimization method in engineering 
and science that enable the individuals of an optimization problem evolves to better solutions j^dj, our proposed 
method can realize a general cloak without following the ideal transformation parameters, but still have a quasi 
perfect performance. Meanwhile, the extreme values exist in the conventional transformation cloak can be avoided. 
When applying the genetic algorithms, the thickness of each layer of the cloak is fixed and the chromosome is a string 
of Os and Is representing a set of constitutive parameters of each layer, {ezi, Mpi, /^^i, ^zMi y^pM^ A'^a/}- The fitness 
of an individual is chosen to be l/Qsca, where Qsca is the total scattering of this individual shown by Eq. ([?]). Making 




FIG. 3; Electric field distributions for a plane wave incident from left to right onto the cloak proposed in Ref. [13 1 



use of the roulette wheel selection, in which individuals with larger fitness have larger chance of going forward to 
the next generation, setting the single point crossover probability and the mutation probability to be 0.6 and 0.05, 
respectively, and ensuring the fittest individual to propagate to the next generation, evolution is carried out and 
optimization is obtained finally. In order to compare the performance of the cloaks created with the TB parameters 
and the optimized parameters, the values of the permittivity and permeability in the optimization are confined to be 
positive values no larger than the maximums of the TB parameters. 

Since in fabrications, it is much easier to create the cloak with fewer layers of matematerials. We consider the 
4-layer cloak as an example. The dimensions of the cloak are Rq ^ \ and R4 = A/2 with a thickness of A/8 for each 
layer. When a TE-polarized plane wave normally incident onto a PEC cylinder with radius A/2 the far-field total 
scattering efficiency Qsca is found to be 2.31. When the PEC cylinder is coated by the 4-layer cloak with the TB full 
parameters, Qsca is reduced to be 0.18. While using our proposed method, we design a 4-layer optimized cloak with 
the parameters shown in Fig. |4](a), (b) and (c), the total scattering efficiency Qsca of this cloak is greatly minimized 
to be 0.0025. 

Fig. [5] shows the calculated electric field distribution when a TE polarized plane wave normally incident from left 
to right onto the optimized cloak. It is shown that in the near region of the cloak, the electric fields stay almost 
unperturbed, thus a so called "quasi-perfect" cylindrical multilayered cloak is obtained. In order to sec how much 
improvement has been achieved through the optimization, the differential scattering efficiency as a function of the 
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FIG. 4: The relative constitutive parameters and the relative impedance for the 4-layer quasi-perfect cloak (solid lines), the 
transformation based cloak (dotted lines) and the ideal linear cloak (dashed lines). 




FIG. 5: Electric field distributions for a plane wave incident from left to right onto a quasi-perfect multilayered cylindrical 
cloak. 

scattering angles is calculated as shown in Fig. [6l We see the scatterings are minimized by about lOdB when the 
TB full parameter 4-layer cloak (dotted line) is used comparing with a bare PEC core (dash dotted line). When the 
4-layer cloak is created with the optimized parameters (solid line), about 20dB of reduction are obtained for all the 
directions comparing with the TB one. Thus a "quasi-perfect" invisibility cloak is achieved although only 4 layers are 
used. 

The reason why quasi-perfect invisibility can still be achieved in a cloak with only a few layers can be physically 
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FIG. 6: Comparison of the differential scattering efficiency of a TE wave incidence for the PEC cylinder with radius A/2 (dotted 
line), the same PEC cylinder coated by the TB full parameter multilayered cloak (dashed line) and by the quasi-perfect cloak 
(solid line). 

explained as follows: The scattering of such kind of multilayered cloak is determined by the recurrence equation Eq. 
[31 When m = 0, Eq. [3] indicates that the nth order scattering is the sum of the direct scattering at Rq denoted by 
roiri and the wave coming out from Rq caused by the multiple reflections and transmissions on the inner boundaries, 
denoted by tion. In order to minimize the total scattering, the parameters of the cloak should be chosen so that 
roin and iion cancel each other. And as denoted by Eq. [31 foi„ are actually determined simultaneously by Vmm km 
and rjm in all the layers. Thus the match of impedances between the conjoined layers does not assure a small total 
scattering, while on the contrary, the mismatch of the impedances, as shown by Fig. [H (d), can be utilized to form 
multiple reflections and transmissions among the inner layers, which eventually produce a transmission to the free 
space being able to destructively interfere with the direct reflection occurs at the outer boundary. In our proposed 
quasi-perfect cloak, the 0th direct reflection coefficient at Rq is tqio = —0.10741 + 0.3096i, while the 0th transmission 
coefRcient from i?o is toio = 0.10739 — 0.3052z, therefore, the 0th scattering coefficient is foio = —0.00002 -I- 0.0044i 
which is very small, due to the destructive interference. This is also similar to the one-dimensional multilayered case, 
where the impedances of each layer are not necessary to be matched in order to get zero reflection. 

It's interesting to sec from Fig.[4]that the optimized parameters (solid lines) are quite different from the parameters 
obtained by the coordinate transformation (dotted lines). Taking the p component as an example, as shown in Fig.[31 
(a), the value of relative permeability in the most inner layer is optimized to be 0.58, instead of a close-to-zero value 
as suggested by the transformation method. We also see that the parameters achieved here are relatively small and 



within the hmit of mctamatcrials, which shows the possibihty of reahzing such a "quasi-perfect" cloak. This is a very 
important contribution to the implementation of the cloak in practice. As we know, for an ideal cylindrical cloak, 
the (p component of the parameters will go infinity near the inner boundary, as shown by the dashed line in Fig. [4] 
(b). Such kind of extreme value near the inner boundary is very difficult to realize. A truncation method at the 



inner boundary can be used to avoid the extreme value of the inner boundary of the cloak, however, scattering will 
be aroused and it has been shown that the performance of such a cloak is sensitive to the perturbations on the inner 



boundary In order to get a performance as good as our proposed cloak, here Qsca is 0.0025, only about A/10^ 
truncation is allowed on the inner boundary, which means huge values of permeability about 5 x 10^/io is needed 
in the direction near the inner boundary, this is a disadvantage of the transformation based cloak in practical 
implementation. 

In conclusion, the analytical model for the multilayered cylindrical cloak has been well established. Based on this 
model, the effects of discretization and simplification, as well as losses, of the transformation based cloak can be 
quantitatively characterized. By utilizing the genetic algorithms, we further show that, it is not the best way to 
manually assign the transformation parameters for the multilayered cloak. A 4-layer "quasi-perfect" cylindrical cloak 
is proposed, whose parameters are relatively small, and possible for realization. The parameters obtained do not 
follow the method of coordinate transformation. Our method was shown to be effective in analyzing a multilayered 
cloak and provides a robust way of designing a practical cloak. 

This work is supported by the ONR under grant No. N00014-06-01-0001, the Department of the Air Force under 
Air Force Contract No. F19628-00-C-0002, and the NSFC under grant Nos. 60531020 and 60801005. 



[1] J. B. Pendry, D. Schurig, and D. R. Smith, Science 312, 1780 (2006). 

[2] A. Greenleaf, M. Lassas, and G. Uhlmann, Physiological Measurement 24, 413 (2003). 

[3] A. AIu and N. Engheta, Phys. Rev. E 72, 016623 (2005). 

[4] A. AIu and N. Engheta, Phys. Rev. Lett. 100, 113901 (2008). 

[5] U. Leonhardt, Science 312, 1777 (2006). 

[6] D. Schurig, J. B. Pendry, and D. R. Smith, Opt. Express 14, 9794 (2006). 

[7] S. A. Cummer, B.-I. Popa, D. Schurig, D. R. Smith, and J. B. Pendry, Phys. Rev. E 74, 036621 (2006). 

[8] F. Zolla, S. Guenneau, A. Nicolet, and J. B. Pendry, Opt. Lett. 32, 1069 (2007). 

[9] H. Chen, B. I. Wu, B. Zhang, and J. A. Kong, Phys. Rev. Lett. 99, 063903 (2007). 



9 

[10] B. Zhang, H. Chen, B. I. Wu, Y. Luo, L. Ran, and J. A. Kong, Phys. Rev. B 76, 121101(R) (2007). 
[11] Z. Ruan, M. Yan, C. W. Neff, and M. Qiu, Phys. Rev. Lett. 99, 113903 (2007). 
[12] B. Zhang, H. Chen, B.-I. Wu, and J. A. Kong, Phys. Rev. Lett. 100, 063904 (2008). 

[13] D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. Pendry, A. F. Starr, and D. R. Smith, Science 314, 977 (2006). 
[14] W. S. Cai, U. K. Chettiar, A. V. Kildishev, and V. M. Shalaev, Nature Photonics 1, 224 (2007). 
[15] M. Yan, Z. Ruan, and M. Qiu, Phys. Rev. Lett. 99, 233901 (2007). 

[16] H. Y. Chen, Z. X. Liang, P. J. Yao, X. Y. Jiang, H. R. Ma, and C. T. Chan, Phys. Rev. B 76, 241104(R) (2007). 
[17] Y. Huang, Y. Feng, and T. Jiang, Opt. Express 15, 11133 (2007). 
[18] J. Roth, and M. J. Dignam, J. Opt. Soc. Am. 63, 308 (1973). 

[19] W. C. Chew, Waves and Fields in inhomogeneous Media, 2nd ed, IEEE Press, 1995. 

[20] D. A. Coley, An introduction to geneitic algorithms for scientists and engineers, World Scientific Publishing Co. Pte. Ltd., 
1999. 



